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ABSTRACT 

The results from a systematic study of eleven pulsar wind nebulae with a torus structure ob- 
served with the Chandra X-ray observatory are presented. A significant observational correlation 
is found between the radius of the tori, r, and the spin-down luminosity of the pulsars, E. A 
logarithmic linear fit between the two parameters yields logr — (0.57 ± 0.22) logE -22.3 ± 8.0 
with a correlation coefficient of 0.82, where the units of r and E are pc and ergs s"""^, respectively. 
The value obtained for the E dependency of r is consistent with a square root law, which is 
theoretically expected. This is the first observational evidence of this dependency, and provides 
a useful tool to estimate the spin-down energies of pulsars without direct detections of pulsation. 
Applications of this dependency to some other samples are also shown. 

Subject headings: pulsars: general — stars: neutron — shock waves — X-rays: stars 



1. Introduction 

Active pulsars eject relativistic pulsar winds 
comprised of relativistic particles and magnetic 
field. Such winds are terminated by a strong shock 
where pressure balance is attained with the ambi- 
ent medium. High energy particles diffusing out 
in the downstream of the shock emit radio to very 
high energy gamma-rays via synchrotron or in- 
verse Compton processes, which are observed as 
pulsar wind nebulae (PWNe). The distance from 
the pulsar to the termination shock, , is expected 
to be 

1/2 

I , (1) 



E 



where E, c, 77, Poxt are the spin-down luminosity of 
the pulsar, the light s peed, the filling facto r, and 
the external pressure (jRees fc GunnI 119741 ). The 
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size of shocks should be rough ly 0.01-0.1 pc with 
typic al parameters of PWNe (jKennel fc Coroniti 
19841 ). In the case of a weakly-magnetized pul- 



sar wind, which is believed to apply to most of 
these systems, particle energy density dominates 
that of magnetic field energy at the shock in which 
the emissivity of synchrotron radiation is relatively 
low. Then, the shock is generally invisible with the 
current instruments exce pt for the brightest ex - 



ample, the Crab Nebula (jWeisskopf et al. I l2000h . 



The equipartition between particle and magnetic 
field energy, at which the Synchrotron emissiv- 
ity is the highest, is reached at a distance of a 
few times of rg, as the po st shock flow decelerates 



(jKennel fc Coromtilll984l ) . Considering that most 
of the pulsar wind energy is blown near the equa- 
torial plane and that Synchrotron cooling becomes 
efficient at the outer region, a torus-like structure 
is expected in this system whose radius is a few 
times of Tg. 

The Chandra X-ray observatory has actually 
revealed the torus (and jet-like) structures from 
several PWNe thanks to the excellent spatial res- 
olution of 0.5 arcsec (c.f., iKargaltsev fc Pavlov 
i2008l ). iNg fc Romalfil (|2004l ) developed a sophis- 



ticated method to fit three-dimensional model to 
the tori, which t hey applied to the Ch andra data 
of several PWNe (|Ng fc Romanill2008al ). providing 
most reliable values of the torus radius. Since the 



1 



termination shock radius is expected to be pro- 
portional to the square root of the spin-down lu- 
minosity, the torus radii may scale in the same 
way as the shock radii, if so, one would expect a 
correlation between the torus radii and the spin- 
down luminosities following the square-root law. 
In this paper, we confirm this relationship for the 
first time, and show its astrophysical use for some 
PWNe. 

2. Samples 

In this paper, we used 11 samples for our study 
with the following cri teria. Chandra observed sev - 
eral tens of PWNe (jKargaltsev fc Pavlov! liooi ). 
Thanks to the excellent spatial resolution of the X- 
ray telescope, equatorial tori and polar jet struc- 
tures are discovered i n mo r e than 10 PWN sys- 
tems. iNg fc Romanil ()2004L l2008al ) measured the 
torus radiu s of 10 PWNe with the developed 
metho d bv iNg fc Romani (2004). Romani et al 



2005 1) uses the same method to B1706 — 44, and 
we added this PWN to our sample. 

Table [T] shows details of our samples. Our 
samples are so young that the pulsar wind is 
still strong and morphological distortion due to 
pulsar motion is relatively small yet; old sys- 
tems sometimes show cometary structure con- 
trolled by this efi'ect. Actually, they are cate- 
gorized to "PWNe with to roidal components" in 
Kargaltsev fc Pavlov! (|2008! ) ex cept for J0537-69 10, 
which has a cometary nebula (jChen et al. Illooi). 
Ng fc Roma^ (|2008a! ) derived the size of its 
torus (farcscc in Table H]) after subtracting the 
diffuse nebula , and w e use the parameters by 
Ng fc Roma^ (|2008ah . The latest distance es- 
timates are used to calculate the physical radius 
of the tori (r in Table[T]) with references in Table[TJ 

3. Results 

In this section, we search for a correlation be- 
tween the torus radius r and spin-down luminosity 
E. Figure [1] represents the plot of r as a func- 
tion of E. Uncertainties of r in the figure are ob- 
tained simply by multiplication of the nebula dis- 
tances and the statistical u ncertainties in appar- 
ent tor us ra dn determined by Ng fc Romani ( 2004 
2008al) and !Romani et al] (|2005l ). It is apparent 
that there is a strong positive correlation between 
log-E and logr, yielding a correlation coefficient 



of 0.82. This power-law like correlation is as ex- 
pected as we review in § 1. On the other hand, it is 
also obvious that there is a non-negligible fluctua- 
tion beyond a simple power-law function owing to 
the statistical uncertainties. These facts suggest 
that r is certainly a function of E but that there 
are other hidden parameters which give a larger 
fluctuation to this relation than the given statisti- 
cal uncertainties. We will discuss possible origin of 
the parameters later. We thus flt the data weight- 
ing them equally with a power-law function, 

logr ^ alogE + (3 , (2) 

where a and /3 are constant values. We obtained 



a = 0.57 ±0.22 
p = -22.3 ±8.0 



(3) 



respectively. The best-fit model is shown in Fig- 
ure H] with a thick solid line. 

In order to check our result from a different 
point view, we calculated the correlation coeffi- 
cient between log£^ and r/E", which should be 
with the best-fit a. Figure 121 shows the relation 
between a and the correlation coefficient. We can 
see that the correlation coefficient becomes close 
to when a is between 0.5 and 0.6. This result 
indicates again that a is around 0.5-0.6. 

The value of a agrees well with the expected 
0.5 (eq.([T])). It also implies that the torus radius 
scales in the same manner as the radius of the ter- 
mination shock. If we fix a to be 0.5, the equation 
becomes 



logr = 0.5 log 19.6 ±0.2 



(4) 



This result is also shown in Figure [T] with a thick 
dashed line. 

There are three samples which are well be- 
low the best-fit lines: Vela, PSR J2229+6114, 
and PSR J1124-5916. The former two are 
known to show clear evidences of interaction 



with ejecta or interstellar medium (jLaMassa et al 



2008; !Kothes et al.!!2001 ).In such a case, as is in- 
dicated by Eq. ^ , the torus radii could be smaller 
confined by higher external pressure compared 
to those without interaction with surrounding 
medium. As for the last one. Park et al. (2004) 
showed that the reverse shock has not yet begun 
to interact with it. However, the latest 510 ksec 
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Chandra observation revealed the almost 3 times 
larger torus compared with that seen in the pre- 
vious 50 ksec Chandra observation which Ng & 
Romani (2008a) analyzed (Park et al. 2007). 

Taking these facts into account, we made an- 
other fit where Vela and PSR J2229-)-6114 are ex- 
cluded and the original value of PSR Jl 124— 5916 
is multiplied by a factor of 3. In this fit, we found 
the better correlation factor of 0.93 and obtained 
the best-fit model of 

logr = (0.49 ± 0.12) \ogE - 19.2 ± 4.7 (5) 

as shown with a thin solid line in Figure [T] 

We also searched for correlation between the 
torus radius and other physical parameters, such 
as age of PWNe, magnetic field, and so on, but 
could not find any significant correlation. 

4. Discussion 

4.1. Termination shocks and tori 

Eq.(l3]) shows that the tori radii (r) show a 
square root dependance on E like the termination 
shock radii (rg): This is the first clue of the termi- 
nation shocks in PWNe. 

Here, we introduce z = r/rs, the ratio between 
the radii of observed torus and the shock. The 
value of z can be common among pulsars in spite 
of different pulsar parameters and environment. 
Assuming eq.([4]), 



77\l/2 



1/2 



1.9 



-0.7 



P.: 



1.6 X 10-9[g cm-is-2] 



1/2 



(6) 



When the temperature, number density, and fill- 
ing factor of external environment of the termi- 
nation shock is 1 keV, 1 cm~^, and 1, {Pext = 
1.6 X 10^^ g cm^^s^^), the ratio of radii of X- 
ray torus and termination shock is almost unity. 
The X-ray emission should come from the out- 
side region of the termin ation shock according to 
Kennel fc Coronitil (1984), which supports our re- 
sult. 

The best-fit a, 0.57, is slightly larger than the 
theoretical value in eq. ([T]) . The data scatter seems 
to be larger at lower values for spin-down energy. 



It could be due to an age effect. The fitting with- 
out Vela and J2229-f6114, which interact with 
ejecta or interstellar medium, shows larger radius 
and smaller dispersion. This result supports our 
scenario. 

4.2. Fluctuation of thermal parameters 

Figure [T] shows that the data points have rather 
large scatter to the best-fit line. We here discuss 
on the possible origin of the fluctuation. Since 
E is derived solely from timing information which 
generally has much higher accuracy than other ob- 
servable quantities, we concentrate on other pa- 
rameters, such as thermal parameters like r, rj and 
Pext, and distance. 

We first estimate how large fluctuation of r is 
required to reproduce the correlation coefficient 
of 0.82 by a simple simulation. A single trial of 
the simulation generates 11 samples from a parent 
population with fluctuated environment, which is 
simulated by log-normal variation in the product 
rjPext with a given standard deviation of ctq- We 
made 1000 trials with a given ctq to calculate an ex- 
pected correlation coefficient < A > and the prob- 
ability for A to be larger than 0.82. Table [5] lists 
the expected correlation coefficient and the prob- 
ability with different cto's- The simulation suggest 
that the observed scatter is explained if ctq ^ 0-6, 
i.e., if fiuctuation in rjPext is a factor of ^ 4. 

One may think that the distance uncertainty 
may be a primary source of this fluctuation. If 
it is the case, we need a factor 2 fluctuation of 
distance since the fluctuation of the torus radii 
is linearly connected to the distance uncertainty, 
which is the square root of rjPext fluctuation. Al- 
though distance to astronomical bodies is not al- 
ways constrained very well, the factor of 2 appears 
to be too large to be attributed to the distance 
uncertainty alone especially for such famous and 
well-studied samples. One of the most famous and 
general measurement method of distance to super- 
nova remnants (SNRs) is to use t he E — D relation 
(e.g.. Case fc Bhattacharvall99"8 ): the distances to 
more than 200 SNRs are estimated using the sur- 
face brightness at 1 GHz and diameter relation, al- 
though there is ^^40% dispersion between distance 
from their method and those from other meth- 
ods. This is because E — D relations can be used 
to estimate properties of ensembles of SNRs, not 



for individual one, as Case fc Bhattacharva ( 19981 ) 
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mentioned. Our relation on the PWN tori can be 
used in similar way to E — D relation to estimate 
their distance. 

The fluctuations of and Poxt result in that 
of Ts, and thus that of r as well. A fluctuation 
of r by a factor of 2 roughly corresponds to that 
of TyPcxt by a factor of 4. This level of fluctua- 
tion could easily occur about Pcxt- The density 
of external gas, which is inside of SNR shells, dif- 
fer by about three orders of ma gnitude from SNR 
to SNR, from ~ 0.1 (SN 1 006; lYa maguchi et alj 
[2OO8.) up to ~ 200 (Gas A; iLazendic et al . 2006,) . 
Thus, we believe that the variation of Pext is the 
primary cause of the scatter. Actually, the torus 
of Vela X, which evidently interacts with ejecta, 
has smaller radius relative to the best fit func- 
tion, which may be due to the larger Poxt effect 
among our samples. Samples without these PWNe 
shows clearly tighter correlation. The fluctuation 
of 77 Poxt is also estimated with eg.® to be 2.5, 
as shown in Table [3l which is much smaller than 
those with all samples. 

The interesting thing is that all of the samples 
with small torii are categorized into combined type 
SNRs, which have radio shell and PWNe. It could 
be due to that SNRs with bright shells have higher 
densities inside the remnants. The detailed sys- 
tematic and observational study of Poxt for indi- 
vidual objects might give a better explanation of 
this fluctuation although it is beyond the scope of 
this paper. 

All we mention above are the "negative" factors 
which tend to wash out the correlation between E 
and r. Nonetheless, we find it with relatively high 
significance, which in turn suggests the robustness 
of this correlation. 

5. Application 

Once the correlation between E and r is es- 
tablished, it will provide a useful tool to estimate 
the spin-down energy of pulsars without direct 
detections of pulsation, although the large error 
range prevent us from precise parameters determi- 
nations. We show some examples of the applica- 
tion. In this section, we cite eq.((3|) although using 
eq.(5) does not alter the results here so much. 



5.1. E determination for GO. 9+0.1 

With eq.lIS]), we can estimate the spin-down lu- 
minosity of PWNe from the radius of tori, without 
information of pulsation. 

Let us consider the example of GO. 9+0.1. 
The SNR GO.9+0.1 has a X-ra y bright PWN 



i n the Galactic center region. iGaensler et al 



(|2001[ ) resolved the PWN with Ghandra. The 
size and flux are 5 arcsec x 8 arcsec and 6.6 x 
10~^^ ergs cm~^s~^ in the 2-10 keV band, respec- 
tively. The torus radius should be about half of 
the longer extension, 4 arcsec, or 0.16 pc with the 
distance of 8.5 kpc. With this value and eq.([3l 
we can estimate the logP of GO. 9+0.1 to be 37.7. 
There is no difference whether we use eq.Q or 
eq.©. 



Possenti et al.l (|2002l ) discovered the empirical 
relation between E and the 2-10 keV luminosity 
Lx (ergs s"!) of PWNe of 



logPx = 1.34 logP - 15.34 



(7) 



altho ugh the dispersion is large (c.f. lKargaltsev fc Pavlov 
I2OO8I ). We can estimate the spin-down luminosity 
of GO.9+0.1 independently from their relation to 
be log E = 37.4, which shows very good agreement 
with our estimate. 

Very recently, iGamilo et al. ( 2009f ) discovered 
coherent pulsation from the central pulsar with 
the period and the period derivative of 52 ms and 
1.5557 X 10"^'^ s s~^, respectively. The resultant 
spin-down energy is log P of 37.6, which also shows 
good agreement with our result. 

5.2. E determination for G328.4+0.2 



Gelfand et al.l (|2007D found an extended struc- 
ture in the PWN G328.4+0.2 with XMM-Newton 
with the size of ^ 1 arcsec, or 0.09 pc at 17 kpc 
(jGelfand et al.ll2007[ ). The coherent pulsation has 
not been detected yet. We can estimate the spin- 
down luminosity in the same way as for GO.9+0.1 
to be log E — 37.3. This is totally consistent wit h 
previous esti mation using 



logP = 37.2 (jGelfand et al 



Possenti et al 
2007h . 



(12002), 



5.3. Distance determination of PSR J1846 0258 

PSR J1846-0258 in SNR Kes 75 is one of the 
mysteriou s PWNe with soft-Gamma-ray repeater 
like flares (jOavriil et allbOOStlKumar fc Safi-Harbl 
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20081 ). The spin-down luminosity is rather large, 
log-B of 36.91. However, the distance of this in- 
teresting source is still unclear: iBecker fc Helfandl 
e stimated that t he dis tance is 19-21 kpc. 



whereas iLeahv fc TianI (|2008l ) suggests that this 
system is much nearer, 5.1-7.5 kpc. 

Eq.(I2]) and the spin-down luminosity indicate 
a torus radius of 0.07 p c, whereas the detected 
torus bv iNg et al.1 (|2008bl ) has the radius of 10 arc- 
sec. It suggests that the expected distance from 
E — r relation is about 1.1 kpc. It is too small 
compared with other distance estimates, or in 
other word, the torus radius is too large. This 
would indicate that either an exceptionally low 
ambient pressure or the pulsar provides addi- 
tional pressure from the inside. It could be 
important information to understand the origin 
of magnetars, which is still hotly debated (e.g., 
Vink &: Kuiped 2006 : Ferrario fc Wickramasinghe 



2006 : iDuncan fc Thom pson' '1993 iGavriil et al 
20081 ) . Systematic study of PWN torii of magne- 
tars should b e done although we have few sample s 
until now (|Rea et al.ll2009l : Ivink fc Bamball2009h . 



5.4. Torus search of DEM L241 

The SNR DEM L241 in the Large Magellanic 
Cloud (LMC) has a compact X-ray source in its 



cente r detected with XMM-Newton (jBamba et al 
20061 ). The flux and photon index are 5.0 x 
10-12 in t;he 2.0-10.0 keV band and 1.57, respec- 
tively. The spin-do wn luminosity expected by 
Possenti et~al1 (|2002f ) is logEdot = 38.4, which is 



one of the largest values among known PWNe. We 
can confirm it when we can detect the torus of 
the PWNe. The size estimated to be 0.39 pc, or 
1.6 arcsec using the distance to the LMC of 50 kpc 
(jPeast 1999,), which can be detectable with excel- 
lent spatial res olution of Chandra , but not with 
XMM-Newton (jBamba et al.ll2006l ). 



6. Summary 

We have made a systematic study of PWN spin 
down luminosity and tori radii using Chandra data 
of 11 samples. It is discovered, for the first time, 
that log r and log E have very strong positive cor- 
relation. The tori could be X-ray emission origi- 
nating from the outside of the termi n ation shocks 
as suggested by iKennel fc Coroniti ( 1984[ ). The 
fluctuation in the correlation is mainly producted 



by variation of the external pressure and distance 
uncertainties. With this correlation, we can esti- 
mate the spin-down luminosity and distance to the 
PWN without information of coherent pulsations. 
This estimation has been applied to GO. 9-1-0.1, 
G328.4-f0.2, PSR 1846-0258, and DEM L241. 
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for Young Scientists. The work of KM is par- 
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Table 1 

Parameters of Tori of PWNe 





radius (rarcscc) 


distance 


radius (r) 


logB 


logT 


References 




[arcscc] 


[kpc] 


[pc] 


[ergs s~^] 


[years] 




PSR J0537-6910 (N157B) 


4.0±0.58 


50 


0.97±0.14 


38.68 


3.70 


(1) (4) 


Crab" 


41.33±0.20 


2 


0.400±0.002 


38.66 


3.09 


(2) (5) 


PSR B0540-69 


2.35±0.06 


50 


0.57±0.01 


38.17 


3.22 


(1) (4) 


PSR J1833-1034 (G21.5-0.9) . . . 


5.7±0.2 


4.8 


0.133±0.004 


37.52 


3.69 


(1) (6) 


PSR J0205+6449 (3C58) 


16.3±0.3 


3.2 


0.253±0.005 


37.43 


3.73 


(1) (7) 


PSR J2229+6114 


9.3±0.4 


0.8 


0.036±0.001 


37.35 


4.02 


(2) (8) 


PSR J1124-5916 


0.9±0.1 


6 


0.026±0.003 


37.07 


3.46 


(1) (9) 


SNR G54.1+0.3 


4.6±0.1 


8 


0.178±0.004 


37.06 


3.46 


(2) (10) 


Vela 


21.25±0.50 


0.29 


0.030±0.001 


36.86 


4.05 


(2) (11) 


PSR B1706-44 


3.30±0.13 


3 


0.048±0.002 


36.53 


4.23 


(3) (12) 


PSR B1800-21 


3.1±0.5 


3.8 


0.057±0.009 


36.34 


3.21 


(1) (13) 



Note.— fl') INg fc Romanil l l2008aD: (2) INg fe Romanil ||200 4|): (3) IRomani et a.1.1 1 12003): (4) IFeasti ||1999 

Trimble fc Woltierl Ill97r); (6) ITian &: Leah-v| 1120081) : (1 \ iRoberts et al.l II1993I) : (8) [Kothes et al.l I I2001D 
Gaensler fc Wallacej 12 003): 
Cordes fc Laziol ||2003) 



(10) IKoo et all I2OO8I) : (11) ICaraveo et al.l l|2001l) : (12) IRomani et al.l ||2005 



(5) 
(9) 
(13) 



^We used parameters for outer torus. 



Table 2 

(To vs. THE expected CORRELATION FACTOR AND THE PROBABILITY OF THE REQUIRED CORRELATION 

FACTOR with eq.dni) 



cro CTo/2 < A> P(A > 0.82) 



0.4 0.2 0.90 0.93 

0.56 0.28 0.82 0.56 

0.6 0.3 0.80 0.49 

0.8 0.4 0.71 0.23 

1.0 0.5 0.62 0.12 



Table 3 

(jo vs. the expected correlation factor and the probability of the required correlation 

factor with eq.® 



cro <To/2 < A> P(A > 0.93) 



0.2 0.1 0.98 1.00 

0.3 0.15 0.95 0.79 

0.36 0.18 0.93 0.56 

0.4 0.2 0.91 0.42 

0.6 0.3 0.83 0.11 

0.8 0.4 0.74 0.04 
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Fig. 1. — Spin-down luminosities vs. radii of tori. 
Triangles and circles are samples with/without in- 
teraction with interstellar medium or ejecta, while 
the box is a new estimation of the radius for 
PSR J1124--5916 (see text). Errors are quoted 
from Table [1] Thick solid and dashed lines repre- 
sent cq.Q and ([4]), whereas the thin solid line is 
for eq.([5]). 




Fig. 2. — Correlation factors with different a. 



